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Abstract — This paper investigates tlie relationship between the 
rank weight distribution of a hnear code and that of its dual code. 
The main result of this paper is that, similar to the Mac Williams 
identity for the Hamming metric, the rank weight distribution of 
any linear code can be expressed as an analytical expression 
of that of its dual code. Remarkably, our new identity has 
a similar form to the MacWilliams identity for the Hamming 
metric. Our new identity provides a significant analytical tool 
to the rank weight distribution analysis of linear codes. We 
use a linear space based approach in the proof for our new 
identity, and adapt this approach to provide an alternative proof 
of the MacWilliams identity for the Hamming metric. Finally, 
we determine the relationship between moments of the rank 
distribution of a linear code and those of its dual code, and 
provide an alternative derivation of the rank weight distribution 
of maximum rank distance codes. 

I. Introduction 

The rank metric has attracted some attention due to its 
relevance to wireless communications [1], [2], public -key 
cryptosystems [3], and storage equipments (see, for example, 
[4]). Due to these applications, there is a steady stream of work 
that focus on general properties of codes with the rank metric 
[4]-[14]. Despite these works, many open problems remain 
for rank metric codes. For example, it is unknown how to 
derive the rank weight distribution for any given linear code 
except when the code is a maximum rank distance (MRD) 
code [5]. Besides the minimum rank distance, the rank weight 
distribution is an important property of any rank metric code, 
and determines its error performance in applications. 

In this paper, we investigate the rank weight properties of 
linear codes. The main result of this paper is that, similar to 
the MacWilliams identity for the Hamming metric, the rank 
weight distribution of any linear code can be expressed as 
an analytical expression of that of its dual code. Our new 
identity is a significant analytical tool for both rank weight 
distribution and hence error performance analysis of linear 
codes. To our best knowledge, no similar result exists in the 
literature. It is also remarkable that our new MacWilliams 
identity for the rank metric has a similar form to that for 
the Hamming metric. Despite the similarity, our new identity 
is proved using a different approach based on linear spaces. 
Using the same approach, we give an alternative proof of the 
MacWilliams identity for the Hamming metric. Based on our 
new identity, we also derive an expression that relates moments 
of the rank distribution of a linear code to those of its dual 
code, and provide an alternative derivation for the rank weight 
distribution of MRD codes. 



The rest of the paper is organized as follows. Section |ll] 
reviews necessary backgrounds in an effort to make this paper 
self-contained. Section IIII-AI introduces the concepts of q- 
product and q-derivative for homogeneous polynomials, and 
investigates their properties. Using these tools. Sections IIII-BI 
and IIII-Cl prove the MacWilliams identity for the rank metric, 
and Section IIII-DI derives the relationship between the mo- 
ments of the rank distribution of a linear code and those of its 
dual code. We also provide an alternative derivation of the rank 
distribution of MRD codes in Section IIII-EI Some examples 
are provided in Section UlI-FI to illustrate our results. Finally, 
Section |IV] adapts the approach in Sections IIII-BI and IIII-CI to 
provide an alternative proof of the MacWilliams identity for 
the Hamming metric. All the proofs have been omitted due to 
limited space, and they will be presented at the conference. 



II. Preliminaries 



A. Rank metric 



Consider an n-dimensional vector x = 
{xn,xi, . . .,Xn-i) e GF((7™)". Assume {ao,ai, . . . ,a,„_i} 
is a basis set of GF(q™) over GF(g), then for 



J =0,1,. 
where a;, 



,n - 



1, Xj can be written as xj 



GF{q) for 



Em — 1 



0,1, 



1. Hence, 



Xj can be expanded to an m-dimensional column vector 
{xQ,j,xi,j, . . . ,Xm~i,j)^ with respect to the basis set 
{ao, ai, . . . , am-i}- Let X be the m x n matrix obtained by 
expanding all the coordinates of x. That is. 



/ 



X = 



xo,i 

XlA 



\ -^rn—l.O *^m — 1,1 



J 



where Xj — J2iLo^ ^ij'^i- The rank norm of the vector 
X (over GF{q)), denoted as rk(x|GF(q)), is defined to be 

def 

the rank of the matrix X over GF((7), i.e., rk(x|GF(g)) = 
rank(X) [5]. In this paper, all the ranks are over the base 
field GF{q) unless otherwise specified. To simplify notations, 
we denote the rank norm of x as rk(x) henceforth. 

The rank norm of x is also the number of coordinates in 
X that are linearly independent over GF(q) [5]. The field 
GF(g™) may be viewed as an m-dimensional vector space 
over GF(q). The coordinates of x thus span a linear subspace 
of GF(g'"), denoted as S(x), and the rank of x is the 
dimension of S(x). 



def 

For all X, y e GF(g™)", it is easily verified that (iR(x, y) = 
rk(x — y) is a metric over GF{q™')", referred to as the rank 
metric henceforth [5]. The minimum rank distance of a code, 
denoted as d^, is simply the minimum rank distance over all 
possible pairs of distinct codewords. 

B. The Singleton bound and MRD codes 

The minimum rank distance of a code can be specifically 
bounded. First, the minimum rank distance of a code 
of length n over GF(g™) is obviously bounded above by 
minjrn, n}. Codes that satisfy d^ = m are studied in [8]. Also, 
it can be shown that < d^ [5], where is the minimum 
Hamming distance of the same code. Due to the Singleton 
bound on the minimum Hamming distance of block codes 
[15], the minimum rank distance of a block code of length n 
in < m) and cardinality M over GF(g™) thus satisfies 

4 < n - log^„. M + 1. (1) 

As in [5], we refer to codes that achieve the equality in Eq. ([U 
as MRD codes. It is also shown that the dual of any MRD 
code is also an MRD code [5]. Clearly MRD codes are the 
counterparts of maximum distance separable (MDS) codes. 

C. Weight enumerator and Hadamard transform 

We restrict our attention to the Hamming metric and the 
rank metric only henceforth in this paper 

Definition 1 (Weight function): Let d be a metric over 
GF(g™)", and define w(v) = (i(0,v) as a weight over 
GF(g™)". Suppose v G GF(g™)" has weight r, then the 
weight function of v is defined as /w(v) = y^x"^'^. 

We shall henceforth denote the Hamming weight function 
and the rank weight function as /„ and /r respectively. Note 
that n is the maximum weight for both /„ and /r. 

Definition 2: Let C be a code of length n over GF{q™'). 
Suppose there are Ai codewords in C with weight i, then the 
weight enumerator of C, denoted as Wc{x,y), is defined as 

n 

W^{x, y) /w(v) = E ^^y'^"''- 

vGC 4=0 

Definition 3 (Hadamard transform [15]): Let C be the 
field of complex numbers. Let a e GF(g™) and let 
{1, ai, . . . , Ofm-i} be a basis set of GF(q™). We thus have 
a = ao + oiai + . . . + a,„_ia„i_i. Finally, let C G C be 

def 

a primitive q-th root of unity. We define — C"°- For 

a mapping / from F to C, the Hadamard transform of /, 
denoted as /, is given by 



(2) 



D. Notations 



In order to simplify notations, we shall occasionally denote 
the vector space GF(g™)" as F. We denote the number of 
vectors of rank m (0 < u < min{m, n}) in GF(g™)" as 
Nu{q"\n). It can be shown that Nu{q"^,n) = ["]a(m, u), 
where a{m,u) is defined as follows: a(m,0) = 1 and 
a{m,u) — n"=o^(9™ ~ 9*) foi" u > 1. The ["] term is the 



Gaussian binomial [16], defined as ["] — a{n,u) / a{u,u). 
Note that ["] is the number of w-dimensional linear subspaces 

of GF(q)". We also define /3(to,0) = 1 and j3{m,u) =^ 
nr=o^ [Tl " > 0, which are used in Section|IinA] These 
terms are closely related to the Gaussian binomial: /3(to, u) — 
\^]f3{u,u) and /3(m + w, m + u) = ["'+"] /3(to, to)/3(w, u). 

III. MacWilliams identity for the rank metric 
A. q-product and q-derivative of homogeneous polynomials 

Definition 4 (q-product): Let a{x,y;m) = 

El=o"*(™)y'^''~' Kx,y\m) = Yfj^Qbj{m)y^x''-^ be 
two homogeneous polynomials in x and y of degrees r and s 
respectively with coefficients ai{m) and bj{m) respectively. 
ai{m) and bj{m) for i,j > in turn are real functions of 
m, and are assumed to be zero unless otherwise specified. 
The g-product of a{x,y;m) and b{x,y;m) is defined to be 

def 

the homogeneous polynomial of degree (r + s) c{x, y; m) — 
a{x, y\ m) * b{x, y; m) = I]uto c„(m)y"x''+''~", with 

Cu{m) = E q"-" ai{m)bu-i{m - i). 

i=0 

For n > the n-th g-power of a{x, y; m) is defined recur- 
sively: a(x, y; m)["l = 1 and a(x, y; m)["l — a{x, y; m)!"^-'^! * 
a(a;, y] m) for n > 1. 

To illustrate the g-product, we provide some examples of the 
q-product. We have x*y — yx, y*x = qyx, yx*x = qyx^, and 
yx*{q"'' — l)y = [q"'' — q)y'^x. Note that x*y ^ y*x. It is easy 
to verify that the g-product is in general non-commutative. 
However, it is commutative for some special cases. 

Lemma 1: Suppose a{x,y;m) = a is a constant indepen- 
dent from m, then a{x,y;m) * b{x,y;m) — b{x,y;m) * 
a(x,y]m) — ab{x,y;m). Also, if deg[c(x, y; m)] — 
deg[a(a;, y; m)], then [a(x, y; m) + c(x, y; m)] * b{x, y; m) = 
a(x, y; m)*b{x, y; m)+c(x, y; m)*b{x, y; m), and b{x, y; m)* 
[a{x, y; m)+c{x, y; m)] = 6(x, y; ■m)*a{x, y; 'm)+b{x, y; m)* 
c{x,y;m). 

def 

The homogeneous polynomials ai{x,y;m) — [x + (q™ — 
and bi{x, y; m) =^ (x — y)['l turn out to be very impor- 
tant to our derivations below. The following lemma provides 
the analytical expressions of ai{x,y\m) and bi{x,y;m). 



def 



Lemma 2: For i > 0, 0"^ 
yW = q°''y'- and x^'l — x^ . Furthermore, 



— . For ^ > 0, we have 



ai{x,y-in) = 



I 

E 

u=0 
I 



a{m, u)y^x 



u^l — u 



bi{x,y;m) = E 



n=0 



(-l)"9-"y"x'- 



(3) 



(4) 



Note that ai{x,y;m) is the rank weight enumerator of 
GF((7™)'. 

Definition 5 (q-transform): We define the q-transform of 
a{x,y;m) = 'Y^\^Qai{rn)y'^x'^~^ as the homogeneous poly- 
nomial d{x,y;m) = X][=o * x^^~^^. 



Definition 6 (q-derivative [17]): For q > 2, the q- 
derivative for x ^ of a real-valued function f{x) is defined 
as 

N def fjqx) ~~ fix) 



{q-l)x 



The g-derivative operator is linear For > 0, we shall denote 
the partial i^-th g-derivative of f{x,y) (with respect to x) as 
f'^'^\x,y). The 0-th g-derivative of f{x,y) is defined to be 
f{x,y) itself. 

Lemma 3: For v <n, the v-t\\ (/-derivative of the function 
a;" is given by j3{n,v)x'^^'^ . Also, the v-\h g-derivative 

of !{x,y) = YJi=ohy'^''^' is given by f^''\x,y) ^ 

Lemma 4 (Leibniz rule): For two homogeneous polynomi- 
als f{x,y) = YJi=ofiy'x''~' and g{x,y) = J2j=o djV^ ^'^^ 
with degrees r and s respectively, the f-th {v > 1) g-derivative 
of their g-product is given by 



{f{x,y)*g{x,y)f''^ 



E 

(=0 



••■ (5) 
Next, we derive the g-derivatives of ai{x,y;m) = [x + 
(q" - and bi{x,y;m) ^ {x ~ y)^'l 

Lemma 5: For h' < I we have 



a] '(x,y;m) 
b'f\x,y;m) 



(3(l,u)ai_^{x,y\m) (6) 
P{l,i^)bi^^{x,y;m). (7) 



B. The dual of a vector 

As an important step toward our main result, we derive 
the rank weight enumerator of (v)^, where v G GF((?™)" 

def 

is an arbitrary vector and (v) = {av : a G GF((j''")}. It is 
remarkable that the rank weight enumerator of (v)^ depends 
on only the rank of v. 

Definition 7: For s > 1 the s-th order B-elementary 
extension of an {n, k) linear code Cq is the {n + s, /c + 

def 

s) linear code defined as Cs = {(cq, . . . , c„+s_i) G 
GF(g™)"+1(co,...,c„_i) - (c„,...,c„+,_i)B G Co}, 
where B is an s x n matrix over GF(g). The 0-th order 
elementary extension of Cq is defined to be Cq itself. 

Lemma 6: Let Cq be an (n, k) linear code over GF{q™) 
with generator matrix Gq and parity-check matrix Hp. The 
s-th order B-elementary extension of Cq is the (n + s,k + 
s) linear code Cs over GF{q'"'') with a generator matrix 



G, = 



Go 







I, 



and a parity-check matrix Hs = 



B 

( Ho I — HqB"^ j. 

Corollary 1: Suppose v = {vq, . . . , G GF{q^)^ has 

rank r > 1. Then C = (v)^ is equivalent to the (n — r)-th 
order elementary extension of an (r, r — 1) linear code with 

It is easy to verify that the (r, r — 1) code with c^r = 2 is 
actually an MRD code as defined in Section Hl-BI 

We hence derive the rank weight enumerator of an (r, r — 
1,2) MRD code. Note that the rank weight distribution of 



MRD codes has been derived in [5]. However, we will use 
our results to give an alternative derivation of the rank weight 
distribution of MRD codes later, and thus we shall not use the 
result in [5] here. 

Lemma 7: For r > 1, suppose = {vq, ■ ■ ■ ,Vr-i) G 
GF(^q™Y has rank r < m. Then the number of vectors in 
Cr = (vr)^ with rank r, denoted as A,._r, depends on only r 
and satisfies Ar^r = Q:(m, r—l)—(7''~^^r-i,r-i- Furthermore, 
the rank weight enumerator of Cr is given by 

Wl{x,y) = g-™ (g" - + (g™ - ~ z/)M} 

The following lemma relates the rank weight enumerator of 
a code to that of any of its s-th order elementary extensions. 

Lemma 8: Let Cq C GF(g'")'' be a linear code with rank 
weight enumerator W^^{x,y), and for s > 0, let {x,y) 
be the rank weight enumerator of its s-th order B-elementary 
extension Cs- Then {x,y) does not depend on B and is 
given by 

{x,y)^ M^c'o {x,y)*[x + {r ~l)vP- (8) 
Combining Corollary [T] Lemma |7] and Lemma |8] the rank 

weight enumerator of (v)^ can be determined at last. 

Proposition 1: For v G GF(g™)" with rank r > 0, the 

rank weight enumerator of £ = (v)^ depends on only r, and 

is given by 

Wl{x,y) = q-'^[[x + {q"^-l)yp + {q^-l)--- 
■■■ (x-y)['-U[a; + (q'"-l)jy]["-^}. (9) 

C. MacWilliams identity 

Using the results shown in Section ITlI-BI we now derive the 
MacWilliams identity for rank metric codes. 

Lemma 9: Suppose that for all A G GF(g™)* and all 
u G i^, we have w(Au) = w(u). For v G GF(g™)", let 
us denote (v)^ as C. Then the Hadamard transform of the 
weight function /w, denoted as /„, satisfies 

W^{x,y) = g-™ [W^{x,y) + (<j" - l)/w(v)] . (10) 
Lemma 10: Suppose v G GF((3'™)" has rank r. Then the 
Hadamard transform of the rank weight function is given by 

A(v) = {x- y)M * [x + (g™ - l)jy]["-'-l . (11) 
Let C be an (n, k) linear code over G¥{q"^), and let 
W2(a;,?/) = X]r=o i'-S weight enumerator 

and W^^{x,y) = X)j=o ^iJ^"''^"""' '■^^ i:?mk weight enu- 
merator of its dual code C^. 

Theorem 1: For any linear code C and its dual code C^, 

W^c- V) = i^^c + - 2; - y) , (12) 

where is the g-transform of W^. Equivalently, 

11 n 



3=0 



(13) 



Also, Bj's can be explicitly expressed in terms of A^'s. 



Corollary 2: We have 



1 



AiPj{i;m,n), 



where 



Pj{i; m, n) 



def 



E 



i 




n — i 


I 




J - K 



(-1) 



(14) 
(15) 



D. Moments of the rank distribution 

Next, we investigate the relationship between moments of 
the rank distribution of a linear code and those of its dual 
code. Our results parallel those in [15, p. 131]. 

First, applying Theorem [T] to C^, we obtain 



i=0 



(16) 

By (/-differentiating Eq. (fT6b ly times with respect to x and 
using the Leibniz rule in Lemma |4] as well as the results in 
Lemma |5] we obtain 

Proposition 2: For < < n. 



n — 1/ 

E 

i=0 



n — I 

V 



A, 



in{k- 



E 



n ^ v 



(17) 



As in [15], we refer to the left hand side of Eq. dTTI i as 
moments of the rank distribution of C. We remark that the 
cases where v = Q and v = n sis trivial. Also, Proposition |2] 
can be simplified if v is less than the minimum distance of 
the dual code. 

Corollary 3: Let be the minimum rank distance of C^. 
If ^ < d', then 



E 



n — I 

V 



(18) 



E. Rank distribution of MRD codes 

The rank distribution of MRD codes was first given in [5]. 
Based on our results in Section lTll-DI we provide an alternative 
derivation of the rank distribution of MRD codes. In this 
subsection, we assume n < m. 

First, we obtain the following results necessary for our 
alternative derivation of the rank distribution. 

Lemma 11: Let {ajj^-^Q and {6i}'^o be two sequences 
of real numbers. Suppose that for < j < ^ we have 



Si=0 [, 



1-3- 



bi. Then for < i < I we have hi 



Based on Corollary [3] and using Lemma [TT] we can derive 
the rank distribution of MRD codes when n < m: 

Proposition 3 (Rank distribution of MRD codes): Let C be 
an {n,k,d^) MRD code over GF(g™) {n < m), and let 
W^{x,y) = X]"=o ^iJ'*^" ' i*-^ rank weight enumerator. 
We then have Aq = 1 and for < i < n — d^. 



We remark that the above rank distribution is consistent with 
that derived by Gabidulin in [5]. 

F. Examples 

In this section, we illustrate Theorem [T] and Proposition |2] 
using some examples. For m > 2, consider the (3, 2) linear 
code Ci over GF(g'") with generator matrix 



Gi = 



1 a 1 
1 a 



where a is a primitive element of GF(g'"). It can be verified 
that the rank weight enumerator of Ci is given by W^^ {x,y) — 



3 



i)r 



Applying Theorem [T] we obtain W^j^ (x, y) = x 
l)y X. We can verify by hand that W^j^(x,y) is indeed the 
rank weight enumerator of C^, which has a generator matrix 
Hi — ( —a 1 ). For Ci, both sides of ( [TtI i are given by 
<?^™, <?" [i] , [I] ), and 1 for 1/ = 0, 1, 2, 3 respectively. 

Note that the results hold when m = 2 < n = 3. 

For TO > 4, let us now consider the (4, 2) code C2 over 
GF(g'") with the following generator matrix 



Go 



v29 



v3« 



C2 is actually a (4,2) MRD code with 4 = 3. Hence, its 
dual code is also a (4,2) MRD code with d^ = 3. 
The rank weight enumerators of both C2 and can be 
readily obtained using Proposition |3] and they are given 
by W^Jx.y) = W^^{x,y) = x^ + [^j (g™ - l)y''x^ + 
|g2m _ 2 - [4] (g™ - 1)} It can be verified that W^^ {x, y) 
and W^i(a;,?/) satisfy Theorem [T] For C2, it can also be 
verified that both sides of ^ are q^™^ [4j^m^ [2]' [f]' and 
1 for 1/ = 0, 1, • • • ,4 respectively. 

Finally, consider the (7,4) code C3 over GF(24) with the 
following generator matrix 



G, 



/ 1 /36 /312 

1 

1 /342 (3^ 

V 1 y 



where /3 is a primitive element of GF(24). Its rank weight enu- 
merator is given by (x, y) = x"^ + IQby^x^ + TibQy^x^ + 
580802/4x3, Theorem [fl indicates that the rank weight enu- 
merator of its dual code is given by W^^ (x, y) ~ x'' + 

3 

465?/3a:4 + 3630y''a;3, which can be verified using exhaustive 
search. It can also be verified that both sides of (fTTT l for C3 
are 21^520192,682752,196416,22416,2772,127, and 1 for 
= 0, 1, • • • ,7 respectively. 



n 

d^ + i 



E(-i: 



df,+i 



IV. MacWilliams identity for the Hamming metric 

In this section, we adapt the approach used in our proof of 
jm{j+i) _ ^\ Theorem [T] to provide an alternative proof of the MacWilliams 
^ identity for the Hamming metric. We first derive the Hamming 
(19) weight enumerator of (v)^, where v is an arbitrary vector. 



Then, using this resuh and properties of the Hadamard trans- 
form, we obtain the MacWilliams identity for the Hamming 
metric. 

Definition 8: For s > 1, the s-th order coordinate extension 
of an (n, k) Unear code Cq is defined as the (n + s, k + s) code 

Cs - {(co, . . .,cn+s-i) e GF(g™)"+^|(co, . . . ,c„_i) £ Co}. 
The 0-th order coordinate extension of Cq is defined as Co 
itself. 

We remark that the s-th order coordinate extension is a special 
case of the s-th order B-elementary extension with B = 0. 

Lemma 12: Let Co be an (n, k) linear code over GF(g'"), 
with a generator matrix Gq and a parity-check matrix 
Ho. Then Cg over GF((7™) has a generator matrix Gg — 

^ ^ ^ and a parity-check matrix Hs = ( Ho | ). 

Corollary 4: Suppose v e GF((7™)" has Hamming weight 
r > 1. Then C = (v)^ is equivalent to the {n — r)-\h order 
coordinate extension of an (r, r — 1, 2) MDS code. 

We hence derive the Hamming weight distribution of an 
(r, r - 1,2) MDS code. Note that [15] gives the Hamming 
weight distribution of all MDS codes. However, that proof 
relies on the MacWilliams identity, and thus may not be used 
here. 

Lemma 13: Suppose = (wq: ■ • • I'l'r-i) G GF(q''")'" has 
Hamming weight r. Then Cr — (vr)^ is an (r, r — 1, 2) MDS 
code whose weight enumerator does not depend on and is 
given by 

The following lemma relates the Hamming weight enumer- 
ator of a code to that of its s-th order coordinate extension. 

Lemma 14: Let Cq C GF((7™)'' be a linear code with 
Hamming weight enumerator W^JyX,y\ and for s > 
let (a;, y) be the weight enumerator of its s-th order 
coordinate extension Cg- Then 

(x, y) = Wl (x, y)-[x + (g" - l)y]^ . (20) 
Combining Corollary HI Lemma [13] and Lemma [141 the 
Hamming weight distribution of C can eventually be deter- 
mined. 

Proposition 4: For v e GF(q™)" with Wh(v) = r, the 
Hamming weight enumerator of £ = (v)^ depends on only 
Wh(v), and is given by 

Wl{x,y) = <z-"{[.T + (g™-l)y]" + (g™-l)... 

••• {x-yr[x + {q"^^l)yT-']. (21) 
Lemma 15: Suppose v e GF(g™)" has Hamming weight 
r. Then the Hadamard transform of the Hamming weight 
function is given by 

/„(v) = {x- yr[x + (g™ - \)yr-\ (22) 
Using Lemma [15] we finally establish the MacWilliams 
identity for the Hamming metric. 

Theorem 2: For any linear code C, we have 

w^c- y) = ^ w^c + ir - x-y)- (23) 



We remark that the MacWilliams identities for the Hamming 
and the rank metrics given in Theorems [2] and [Tj respectively 
have exactly the same form except for the g-transform in 
Eq. ( [T2] |. Note that Theorem [2] is precisely the MacWilliams 
identity for the Hamming metric given by Theorem 13 in 
[15, Chap. 5], although our proof is different from that in 
[15, Chap. 5]. Finally, we remark that Theorem 13 in [15, 
Chap. 5] is a special case of the MacWilliams Theorem 
for complete weight enumerators (see Theorem 10 in [15, 
Chap. 5]). For the rank metric, it is not clear how we can 
adapt the concept of complete weight enumerator to give a 
proof of the MacWilliams identity. 
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